Current-induced transverse spin wave instability in a thin nanomagnet 
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We show that an unpolarized electric current incident perpendicular to the plane of a thin fer- 
romagnet can excite a spin-wave instability transverse to the current direction if source and drain 
contacts are not symmetric. The instability, which is driven by the current-induced "spin-transfer 
torque", exists for one current direction only. 

PACS numbers: 75.75.-(-a, 75.40.Gb, 85.75.-d 



Ferromagnets serve as spin filters for an electrical 
current passing through the magnet: the spin of the 
electrons that are transmitted through a ferromagnet 
becomes partially polarized parallel or antiparallel to 
the direction of the magnetization whereas spin cur- 
rent perpendicular to the magnetization direction is ab- 
sorbed. Spin filtering is the root cause for the "spin- 
transfer torque" , the phenomenon that a polarized cur- 
rent impinging on a ferromagnet affects its magnetiza- 
tion direction Q, The source of the spin po- 
larized current can either be a different ferromagnet, 
or, for a thick magnet, a region of the same ferromag- 
net upstream or downstream in the current flow. The 
"spin-transfer torque" gives rise to magnetization rever- 
sal in ferromagnet-normal-metal-ferromagnet trilayers 
P, |3| , which has been observed experimentally by sev- 
eral groups 0, 0, IE 0i Hj IS ^1 ■ Dynamic manifestations 
of the spin-transfer torque include domain wall motion in 
bulk ferromagnets [ij [ij, [ij, [iJI and the excitation of 
spin waves by polarized currents in ferroma gne tic mul- 
tilayers or wires |llllillli,lllllillil,ll3. InaU 
these manifestations, the current-induced spin torque can 
be distinguished from effects arising from the current in- 
duced magnetic field, the main difference being that spin- 
transfer torque effects depend on the current direction, 
whereas magnetic field induced effects do not. 

In this letter, we show that an unpolarized current can 
also exert a spin-transfer torque on a ferromagnet, even 
if the magnet is so thin that its magnetization direction 
does not change along the current fiow: Although an un- 
polarized current cannot exert a spin-transfer torque that 
changes the over-all magnetization direction, it can cre- 
ate a transverse spin wave instability for sufficiently high 
current densities if the source and drain contacts to the 
ferromagnet are not symmetric. This spin wave insta- 
bility can be identified unambiguously as a spin-torque 
effect because of its dependence on current direction: 
the spin- wave instability is present for one current direc- 
tion and absent for the other. The spin-wave instability 
should lead to a non-hysteretic feature in the current- 
voltage characteristic of the ferromagnetic film that ex- 
ists for one current direction only. In thick ferromag- 
nets, such features have been observed in recent experi- 
ments 0, ll^ . The necessary criterion for the spin- wave 



FIG. 1: Left; Schematic picture of a thin ferromagnetic layer 
(F) in a "nanopillar" geometry. The ferromagnet is connected 
to source and drain reservoirs through normal metal leads (N) . 
Current flow through the ferromagnet leads to a spin accu- 
mulation in the normal leads, as shown by the dotted arrows. 
Right: Cartoon of spin wave generation. Spin diffusion of 
reflected electrons (diffusion paths shown full and dashed) re- 
sults in a spin torque on the ferromagnet's magnetization that 
enhances the spin wave amplitude; the transmitted electrons 
(not shown) damp the spin wave. 



instability, asymmetric contacts to source and drain, is 
gcnerically fulfilled in experiments on nanoscale magnets 
5j. The issue of current-induced spin- wave excitation 
has significant practical relevance for devices based on the 
spin-torque effect in ferromagnetic multilayers. Whereas, 
experimentally, the presence of dynamical phenomena in 
these systems is well established, the precise nature of the 
excitations is not well known. A theoretical understand- 
ing of current-induced spin wave excitation in multilayer 
structures will shed light on the possibly useful applica- 
tion of the current-induced dynamical excitations, e.g., 
as GHz resonators, as well as on ways to avoid spin- wave 
excitation in devices that are designed to exhibit mag- 
netization reversal only. For ferromagnet-normal-metal- 
ferromagnet trilayers, a comparison between theory and 
experiment is obstructed by the interplay of spin wave 
excitations and static magnetization reversal. Our find- 
ing that dynamical phenomena exist already in a single 
ferromagnetic layer allows for a study of spin-wave exci- 
tations in a much simpler geometry in which dynamical 
and static phenomena are well separated. 

The geometry of the system under consideration is 
shown in Fig. ^ ^ nanomagnet of thickness d, small 
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enough to be considered as single domain, is connected 
to source and drain reservoirs via diffusive normal-metal 
leads of lengths L_ and L^, respectively. The leads and 
the ferromagnet have width W ^ d. In the absence of 
an electrical current through the system, the ferromagnet 
has a uniform magnetization, with a direction determined 
by anisotropics and an external magnetic field. A current 
j flows perpendicular to the ferromagnet. (Note that the 
electrical current j points opposite to the electron flow.) 

For a qualitative explanation of the mechanism of the 
spin-wave instability, we note that the passage of an elec- 
tric current through a single ferromagnetic layer creates 
spin accumulations of opposite signs on both sides of the 
ferromagnet, see Fig. ^ where the sign of the spin ac- 
cumulation depends on the current direction and on the 
spin- filtering properties of the ferromagnet. To first or- 
der in the spin- wave amplitude, the spin accumulation 
in the normal metal leads is not affected by the possi- 
ble presence of spin waves in the ferromagnet, as long as 
their wavelengths are much smaller than the spin diffu- 
sion length in the normal metal. Each spin accumulation 
exerts a spin-transfer torque on the magnetization; the 
direction of the torque is to align the magnetization in 
the ferromagnet with the direction of the accumulated 
spins in the leads. Depending on the sign of the spin 
accumulation, such a torque either damps or enhances 
non-uniform spin wave excitations in the ferromagnet. 
The magnitudes of the spin accumulations on each side 
of the ferromagnet are, generically, different, since they 
depend on the spin diffusion length, scattering proper- 
ties of normal-metal-ferromagnet interface, distance to 
reservoirs, etc. Therefore, an unpolarized current re- 
sults, generically, in a net torque that either suppresses 
or enhances non-uniform spin waves in the ferromagnet, 
depending on current direction. The instability occurs 
when the current-induced enhancement of a spin wave 
amplitude overcomes the intrinsic spin-wave damping. 

For a quantitative description of the spin- wave instabil- 
ity, we choose the x axis along the leads, with the origin 
such that the normal-metal-ferromagnet interfaces are 
at X = ±d/2, see Fig. ^ We assume that the ferro- 
magnet is so thin that its magnetization direction varies 
with respect to the transverse coordinates y and z only, 
but not with respect to x. In the diffusive normal-metal 
leads, the electron distribution is described by potentials 
lie{^,t) for the electron density and /Xg(r,t) for the elec- 
tron spin |18>. .19.„ 20]. The current is separated into the 
charge current density ja, a = x,y,z, and the spin cur- 
rent density ]sa- Since electrons adjust to the changing 
magnetization on a time scale much faster than that of 
the magnetization dynamics, current and potential are 
related by means of the time-independent diffusion equa- 
tion in the normal metal leads, 

V^/Lte = 0, Jq = (cr/e)VQ^e, a = x,y,z, 



supplemented by boundary conditions for the source 
and drain reservoirs fj,e{—L-) = —eV, iJLe{L+) = 0, 
^ig{±L±) — 0. Here a is the conductivity of the nor- 
mal metal leads and Zgf the spin-diffusion length. At 
the normal-metal-ferromagnet interfaces, the charge and 
spin current jx and }sx in the normal metal perpendicu- 
lar to the interface are related to the potential drop A/i 
over the interface as [ifll EH] 

{}sx)± = {h/2e^)'Reg^[ C^^fJ-s x m±hdtm) x m 

+ (n/2e^)Im.g|| (2A/x^ x m ± hdtm) , 
(j.x)ll = -(.9TT+3ii)(^/2e')m-A/.i, (2) 

- (.9TT - .9ii)(^/2e')AMe, 
jx = (giT +3ii)(AMe/e) + (.9TT - 5u)m • (A/^Je). 

Here m(r, t) is the unit vector pointing in the direction of 
the magnetization of the ferromagnet, (jsx)± and (jsx)|| 
are the z-components of the spin current perpendicu- 
lar and parallel to m, respectively, and are in- 
terface conductivities for spins aligned parallel and an- 
tiparallel to m, whereas gn is the "mixing conductivity" 
[laEi, A/i = fi{±d/2 + 0)~^{±d/2-0) is the potential 
drop over the interface. (Wc assume identical conduc- 
tivities and spin-flip lengths in the leads, and identical 
scattering properties of the ferromagnet-normal-metal 
interfaces, although our results are readily generalized 
to unequal values of cr, 1st and gap-) At the normal- 
metal-ferromagnet interface, (js2;)|| and jx are continu- 
ous, whereas {jsx)i. = in the ferromagnet. We assume 
that the ferromagnet is so thin that all potential drops 
occur at the interfaces, so that we can neglect the x- 
dependence of the potentials in the ferromagnet. Then 
Eqs. and Q fully determine the potentials and cur- 
rents in the normal metal and the ferromagnet, as a func- 
tion of m. 

For a thin ferromagnet, variations of the magnetiza- 
tion direction m along the x axis will have a large energy 
cost and, hence, a large threshold for excitation. There- 
fore, we consider transverse modulation of m only and 
take m independent of x. The dynamics of m{y,z,t) 
is determined by the Landau-Lifschitz-Gilbcrt equation 
(without an applied magnetic field [2, 13], 

dttn ~ am x dtm + JjAIV'^m x m 

- {j/M){Kimiei + K2m2e2) x m (3) 
+ {-f/Md)%xi-d/2 - 0) ~3sxid/2 + 0)], 

where a and J are the bulk Gilbert damping coeffi- 
cient and spin stiffness (exchange constant) respectively, 
7 = fisg/fi is the gyromagnetic ratio, M is the magne- 
tization per unit volume, and Ki and K2 are anisotropy 
constants along principal directions ei and 62, respec- 
tively, obtained by expanding the magnet's free energy 
around a preferred axis 63. The magnetization satisfies 
the boundary condition (n- V)m = 0, where n is the nor- 
mal to the ferromagnet's surface. We neglect the effect of 
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the current-induced magnetic field on the magnetization 
dynamics, which is allowed if the width W is sufficiently 
small, <tC 1 /im for typical experimental parameters. The 
parameters in the Landau-Lifschitz-Gilbert equation de- 
fine a length scale l/q{ and current scale jf, 



The quantities 1/qi and hjf/e are proportional to the 
width and energy of a domain wall, respectively. Order- 
of-magnitude estimates of the various parameters in- 
volved here are d ^ 10 nm, W ^ 10^ nm, /sf(Cu)~ 10^ 
nm Jllll; Img^x < Reg^ ~ 5tT ~ 9U ~ 10^^ f^-^m^^ 
Ullia for Co/Cu and Fe/Cr interfaces; cr/4f (Cu)- lO^s 
51~^m~^, and typically qi ^ 10^^ nm~^, jf ^ 10^ A/cm^ 
for Co, Fe or Ni (2^]. 

We first solve these equations for a uniform magnetiza- 
tion, m independent of the transverse coordinates y and 
z. The spin accumulation /x^ in the normal metal leads 
close to the normal-nietal-ferromagnet interface is 

Ms(±c;/2) = T{ejx/gm)tanh{L±/lsf)m, (5) 

where jx is the charge current density and 



9m = 



(g/^sf)(gTT + + 2.9TT.9iiE±tanh(L±/?sf) 
.9TT - 9 II 



The spin accumulation is shown schematically in Fig. ^ 
In the case of uniform magnetization, the spin accu- 
mulation fig is always parallel to m, and no current- 
induced torque is applied to the magnetization. The sit- 
uation changes if m varies in the transverse direction. 
In this case, transverse diffusion of spin in the normal- 
metal leads gives rise to an angle between fi^ and m, and, 
hence, to a current-induced torque. In order to study this 
scenario in detail, we analyze Eqs. Q-Q for a small 
deviation of m from the equilibrium direction 63. The 
result can be represented in terms of an equation of mo- 
tion for (5m = m — 63. We assume a rectangular cross 
section of dimensions Wy and Wz in the y and z direc- 
tions and perform a Fourier transform with respect to the 
transverse coordinates y and z. The allowed wavevec- 
tors are qy = irny/Wy, qz — nUz/Wz, where Uy and 
Uz are non-negative integers. Representing through 
™±(<l) = Smi{q) ± ii5m2(q), one finds that, to first or- 
der in m-t, the equations of motion of different Fourier 
modes separate. 



a ,i \ n ra K2 - Ki hjijq^ + qf) 

— ± — Al0tm± — —m± 

77/ 2 2eqf 



2ed 



-m±, 



(6) 



Here a and 7 are renormalized Gilbert damping param- 



eter and gyromagnetic ratio, 

7 7^2Mde2 Z-^G±{q)+gu' 



9nG±iq) 



2Mde^ ^G±{q)+gu' 



(7a) 



(7b) 



whereas the dimensionless numbers Si and S2 set the 
magnitude of the current-induced torque, 

S. = ^Rej:^^^^^l^M^^ (8a) 



gm^: 



sf 



G±(0) gn + G±{q) ' 



5. = ^ImV-^^±M-^, (8b) 
gmkf ^G±(0) gn + G±{q) 



where gm was defined below Eq. ((SJ and 



G±iq) 



17/ +q^ COth[L±Jl, 



In the limit q ^ 0, Eqs. (TJ coincide with the enhanced 
Gilbert damping and gyromagnetic ratio reported by 
Tserkovnyak et al. [2l| . 

In the absence of a current, any spatial modulation of 
the magnetization is damped. It is the existence of the 
source terms Si and S2 in Eq. (jSJ that leads to a spin- 
wave instability at sufficiently large current density j. 
Note that the source terms exist only if the normal leads 
are asymmetric (different lengths, or different conductiv- 
ities, spin- flip lengths or interface conductivities), and if 
the ferromagnct is a spin filter, g^^ ^ gn. The source 
term ^2 gives rise to a small change of the ferromagnetic 
resonance frequency, whereas Si increases or decreases 
the amplitude of the spin wave, depending on the current 
direction. An instability occurs if the current-induced 
enhancement of the spin-wave amplitude overcomes the 
damping, i.e., if 



V "7 / 



(g2 + qf)d 



9f 



Jf- 



(9) 



Using the order-of-magnitude estimates listed below 
Eq. Q , the instability criterion simplifies considerably if 
we set Im(7|| = 0, take the limits d and a/lgf ^ 
Re(7||, and consider the case L_ ^ Isi ^ of maxi- 
mally asymmetric contacts. 



jf 



> 



qf + q^ 



Mgfe2l-(l + (;2/2^) 



-1/2- 



(10) 



The r.h.s. of Eq. (|10|l is shown schematically in Fig. |2] 
For large q S> the spin- wave instability is domi- 

nated by the stiffness of the spin wave, which leads to a 
critical current density oc q^. For small q <C 1/^sf, the 
magnetization accumulation in the normal metal leads 
tends to remain locally parallel to the magnetization, so 
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FIG. 2: Schematic picture of the current density je required 
to excite a spin-wave at wavevector q. The spin-wave insta- 
bility occurs at the wavevector q for which j is minimal. 

that the current-induced torque is strongly reduced and 
the critical current density for excitations at small q is in- 
creased oc . For a thick ferromagnet, the effect of bulk 
Gilbert damping becomes dominant, causing the critical 
current density to increase linearly with the thickness d. 

In the experimentally relevant parameter regime 
1/gf ^ ^sf ^ f/Re^Ti' critical current density is 



7gmgf 

Me2 



1/3 



(11) 



For a Cu/Co/Cu geometry, we estimate a ^ (0.1nm)/d 
and ti^jgraqi/Me"^ ^ 0.2, as long as dqi < 1. We then 
find jc ~ lO'' A/cm^, which is of the same order of magni- 
tude as the critical current for magnetization reversal in a 
ferromagnet-normal-metal-ferromagnet trilayerj^. The 
estimate (lllll is valid as long as the width W ^ l/fZc- 
For small W, the instability occurs at the lowest possible 
wavenumber, qc = tt/W. 

In conclusion, we have shown that an unpolarized elec- 
tric current passing through a thin ferromagnetic layer 
can create a spin wave instability with wavevector trans- 
verse to the current. The instability occurs at one current 
direction only and depends on the asymmetry of normal 
metal contacts. The mechanism for the instability is the 
same as the one believed to cause magnetization reversal 
in ferromagnetic multilayers, i.e., a spin-transfer torque 
arising from spin accumulation in the normal metal con- 
tacts perpendicular to the magnetization direction m and 
spin-filtering at the normal-metal-fcrromagnct interfaces 

Current-induced excitations that occur for one current 
direction only have been observed in single ferromagnetic 
layers |j,llall3- For ferromagnets much thicker than the 
inverse wavevector q~^ ~ 3 x lO^nm for transverse spin 
waves, cf. Eq. Hll|l . spin wave excitations with wavevec- 
tor along the current flow may have a lower threshold 
current than transverse excitations 0, 0| . Such thick 
ferromagnetic layers are common in the point-contact ge- 
ometry 0, 0, 01 . It is_foi_the ^in layers used in the 
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nanopillar geometry |5|, [3, la IE |23j that we believe the 
transverse instability considered here is most significant. 



